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In double parton interactions, the two hard proceses are correlated via double parton densities. 
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cross section. In particular correlations between the spins of the interacting quarks can induce 
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correlations following from positivity of the double parton densities. 
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1. Introduction 

The increase in the density of partons with the energy of a hadron collider increases the amount 
of interactions where more than one parton in each hadron participates in hard interactions. Multi- 
parton interactions are therefore expected to play a larger role at the LHC than ever before. They 
contain information about the proton structure not available from single-parton interactions, and 
they can be important backgrounds to other signal processes, such as Higgs production [1]. Multi- 
parton interactions have an experimental history stretching from the ISR [2] to the LHC [3], and 
have long been modeled in event generators [4]. However, the treatment of multi -parton interac- 
tions from the point of perturbative QCD has only recently started to accelerate. 

The hard interactions are connected via parton distributions and the partons originating from 
the same hadron can be correlated. This leads to interferences and correlation effects not present in 
single-parton scatterings. These effects are often neglected in studies of multiple hard scatterings, 
but the validity of such approaches still has to be investigated. A suitable scene for such a study 
is set by the double Drell-Yan process [5, 6, 7], where two quark/anti-quark pairs annihilate into 
two vector bosons (y*, Z, W ± ) which decay leptonically. This process has the advantage of being 
theoretically clean and well understood in the single-parton scattering case. We calculate the dif- 
ferential cross section, taking initial parton correlations into account, and derive constraints on the 
size of the spin correlations following from positivity of double parton densities. 

2. Double Parton Interactions 

When two partons in a proton interact, it is only the sum of momenta and quantum numbers 
which have to be the same in the amplitude and its conjugate. This allows for a momentum differ- 
ence, r, between a parton in the amplitude and its partner in the conjugate amplitude, see figure 1. 
This difference has to be balanced by the parton in the other interaction. Similarly the color and 
flavor of the partons can be matched inside each hard collision or as an interference effect between 
them. There can even be fermion number interference between quarks and anti-quarks but we do 
not include that in our calculations. Further the spin of the interacting partons can be correlated, 
much in the same way as in single-parton collisions with polarized protons [8]. The correlations 
between the initial state partons are described by double parton distributions, DPDs. 

2.1 Double Parton Distributions 

In transverse position space the DPDs depend on the collinear momentum fractions x\ and xj 
of the two partons taking part in hard collisions, on their differences in transverse positions between 
the amplitude and its conjugate z\ and z 2 and the transverse distance between the hard vertices y. zi 
(Z2) is the Fourier conjugate of the average transverse momentum ki (I42) and y of the momentum 
difference r, see figure 1. The DPDs describing the different quark polarizations are labeled by q 
for unpolarized, Aq for longitudinally polarized and 8q for transversely polarized quarks [9]. 

For unpolarized and longitudinally polarized quarks the possible combinations are 

Fqq = fqq(xi,X2,Z U Z 2 ,y), F qAq = g q & q (x\ ,X 2 , Z\ , Z 2 , y) , 
FhqAq = fAqAq (*1 , *2, z l j *2) y) j F&qq = gAqq (*1 ; *2 , z l j z 2 , J) , 
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Figure 1: The double Drell-Yan process where two quarks in the right moving proton interact with two 
anti-quarks from the left moving proton. q\ is the momentum of the vector boson from interaction i. k\ and 
£2 (k\ and k.2) are average momenta carried by the partons from the right (left) moving proton taking part in 
hard interaction 1 and 2 respectively and r the momentum difference. 

where f's are scalar- and g's are pseudo scalar-functions. Distributions with transversely polarized 
quarks in one of the two hard interactions carry an open index corresponding to the transverse spin 
vector, and have to be expanded in a basis spanning the transverse plane 

F L, 8q = M {y'fMS, +f g^sq) , P l q8q = M {ff q8q +y i g q8q ) (2.2) 

and similarly with the subscripts interchanged. M is the proton mass and f = y^e'\ is a transverse 
vector orthogonal to y l . Transversely polarized quarks in both interactions give two transverse 
indices and we need a tensor basis 

= §iJ f^ q + (2yV - y 2 8 ij ) M*f 8q8q + ( y y + fyJ ) Mh s 5q5q 
+ (y i y j -y i y j )MY 8qSq . 

Taking color interference into account gives one color singlet and one color interference distribution 
for each flavor combination and thus doubles the number of DPDs. When the flavors of the quarks 
are different there are flavor square and interference distributions. The DPDs for the left moving 
proton will be denoted by a bar, f qq , which should not be confused with the bar appearing over 
subscripts (indicating anti-particles). 

3. Cross Section 

The partonic cross sections are calculated to leading order for unpolarized, longitudinally and 
transversely polarized quarks and the results are combined with the DPDs describing the corre- 
sponding parton densities. 

The cross section formula, assuming transverse momentum dependent factorization, derived 
in [9], can be expressed as 

YlUdxidxid^idcii c qi fy 3 ^ aia % 4 da, da 2 J (i%y e J {i%y e (31) 

x / d y[F a + 

^0104^2302 ^0302^0104 F a3Cl4 F CIia2 ] + {flavor interference} . 
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Figure 2: Reference frame for each interaction is the rest frame of the vector boson. 0, is the polar and <p, 
the azimuthal angle of the lepton. k (If) is the three momentum of the lepton (anti-lepton) from interaction i. 

C equals 2 when the final states of the two hard interactions are identical and 1 otherwise. In 
the first sum q\ (#2) labels the flavor of the quark and 173 (^4) of the anti-quark entering the first 
(second) hard scattering, a,- and a,- label the different combinations of unpolarized (qi) and lon- 
gitudinally (A<7,-) and transversely (8qi) polarized quarks and anti-quarks. The color square and 
interference distributions enter the cross section with equal weights and we can therefore keep 
color square/interference labels implicit. ddi/dQ-i is the partonic cross section differential in the 
direction of the outgoing lepton. 

To describe the final state we use angles defined in the rest frames of the vector bosons, with 
z-axes defined as the direction bisecting the angle between p and — p (p and p are three momenta of 
the protons), see figure 2. The x-axis is an arbitrary transverse direction which we for definiteness 
choose to point towards the center of the LHC ring. The differences between the axes defined for 
the two vector bosons are of order qt/Qi and can be neglected. 

We will display the final results only for the cross sections integrated over transverse boson 
momenta explicitly, not including flavor interference, and refer to [10] for the cross sections de- 
pending on transverse boson momenta. 

Integrating over the transverse momenta of the bosons yields collinear double parton distri- 
butions [5]. After integration the DPDs depend on x\, x 2 and only one transverse vector, y, and 
therefore no pseudo-scalar functions can contribute. Due to time reversal symmetry, functions 
with one longitudinal and one transversely polarized quark vanish, for example f^Sq = 0. We are 
then left with six collinear double parton distributions for each combination of quark flavors. For 
unpolarized and longitudinally polarized quarks the cross section is 



1 



£ [^(1+ cos 2 00 +4^008 0!] [^ 4 (l+cos 2 2 )+^ 4 cos0 2 ] 



dxjXjdQj 



C 




(3.2) 



+ [K qi Aq 3 ( 1 + cos 2 0i ) + K' qiAii3 cos 0i] [K qiA g 4 ( 1 + cos 2 2 ) + K' qiAq : 4 cos 2 ] 
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K q q, K' qq , K q hq and K' qA = are Q 2 dependent coupling factors for different polarizations, denned in 
[10]. £2,'s are solid angles and 0;'s polar angles of the leptons. 'perm' stands for permutations 
of the quark/anti-quark subscripts in the DPDs. For the coupling factors are polarization in- 
dependent and including partonic spin correlations therefore only changes the rate, while for the 
neutral current the longitudinal polarization changes both the rate and the angular distribution of 
the process. 

The part of the cross section with one interaction involving transversely polarized quarks van- 
ish upon integration over the transverse boson momenta while the cross section with both interac- 
tions containing transversely polarized quarks is non-zero for the neutral current, 



do® 



1 



dxjdxidCli C 



9112 



£ sin 2 0! sin 2 2 / d 2 yl {K 8qi8qi K 8q2 s q2 -K' 8qi8q K' 8q28q2 )co^2{(pi - 92) 



( K S qi S qi K Sg 2 8g 2 - K Sc lt S qi Ksq 2S q 2 )sm2((p l - (fr) 2 [fSq^qJSq^ + PCTm. 



(3.3) 



KsqSq an d K' SqS - are Q 2 dependent coupling factors for two transversely polarized quarks and the 
9,'s are the azimuthal angles of the final state leptons. For W bosons the coupling factors are zero 
for transversely polarized quarks, since the W only couples to left handed quarks. The cross section 
for the neutral current depend on the azimuthal angle between the two outgoing leptons (<pi — (fc). 
This dependence originates in the correlations between the spin of the initial state quarks. 



4. Positivity Bounds 

The double parton distributions of different polarizations can be organized in a positive semi- 
definite spin density matrix. The positivity can then be used to derive upper limits on the polarized 
distributions, similarly to what was done for generalized parton distributions in [1 1] and for trans- 
verse momentum dependent distributions in [12]. 

The projection operators T ++ (r ) project out quarks with positive (negative) helicities while 

r_ + and r + _ give helicity interferences. These can be expressed in terms of operators projecting 
out un- (T q ), longitudinal- (Thq) and transversely- (T 8q ) polarized quarks 



-(i + 75 ) = r,+r Afy , 
-(i-7 5 ) = r,-r Af/ , 



10 



+1 
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2 

+1 



-(1 + 75) 



r 1 



(4.1) 



2 y ~ ,J/ ~ 1 ~ ^Hi ~t 2 v ~ ' " ' ' ' Sq 

The DPDs can then be organized as a matrix in the light-cone helicity basis 

(fqq+fAqAq -fyWfSqq ~i\y\Mf q8q 2y 2 M 2 

i\y\ M fSqq fqq-fhqhq V&q&q ~i\y\ M fqSq 
i\y\Mfq Sq 2f SqS q fqq-fhqhq ~i\y\Mf 8qq 

\2y 2 M 2 f^ qSq i\y\Mfq S q i\y\Mf Sqq f qq + fhqhq) 



Sq' 



(4.2) 



where the rows (columns) correspond to ++, — h, H — , helicities of the two quarks in the am- 
plitude (conjugate amplitude). The non-interference parton distributions can be interpreted as the 
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probability of rinding quarks having specific helicities inside a proton. The helicity matrix is there- 
fore positive semi-definite and calculating the eigenvalues we obtain necessary and sufficient posi- 
tivity conditions 

f qq > \f SqSq ~y 2 M% Sq \ (4.3) 

{f qq ^{f SqSq -y 2 M 2 fl qSq )) 2 -[h qM T(fs qSq WM (4.4) 
This implies the somewhat more transparent but weaker inequalities 

f qq + /a,a ? > 2y 2 M 2 f SqSq f qq - f AqAq > 2f SqSq . (4.5) 

5. Summary 

Spin correlations as well as color and flavor interference result in a large number of double 
parton distributions. Many of these distributions contribute to similar terms in the double Drell- 
Yan cross section. Longitudinal polarization changes the overall rate for the charged bosons, while 
for the neutral bosons the angular distribution is also affected. The spin vectors of transversely 
polarized quarks lead to a dependence on the azimuthal angle between the two outgoing leptons. 
Similar features as those appearing in the double Drell-Yan cross section are expected to be present 
also in other types of processes, such as double dijet production, with larger cross sections but a 
dramatic increase of complexity due to their color structure. 

We used the probability interpretation of double parton distributions to derive constraints on 
the polarized parton distributions. Similar constraints have proven very useful for single parton 
distributions. 
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